Temperature dependence of surface magnetization in local-moment systems 
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We present a theory to study the temperature-dependent behavior of surface states in a ferromag- 
netic semi-infinite crystal. Our approach is based on the single-site approximation for the s-f model. 
The effect of the semi-infinite nature of the crystal is taken into account by a localized perturbation 
method. Using the mean-field theory for the layer-dependent magnetization, the local density of 
states and the electron-spin polarization are investigated at difi'erent temperatures for ordinary and 
surface transition cases. The results show that the surface magnetic properties may differ strongly 
from those in the bulk and the coupling constant of atoms plays a decisive role in the degree of 
spin polarization. In particular, for the case in which the exchange coupling constant on the surface 
and between atoms in the first and second layer is higher than the corresponding in the bulk, an 
enhancement of surface Curie temperature and hence the spin polarization can be obtained. 



I. INTRODUCTION 

Knowledge of the spin-dependent electronic structure 
at surfaces and interfaces plays an increasingly impor- 
tant role when assessing possible use of novel magnetic 
materials for spintronic applications^. A surface breaks 
the translation symmetry of a system and changes lo- 
cal quantities such as local density of states (LDOS) 
which is significantly different from the bulk density of 
states (DOS). Experimental study on magnetic surfaces 
has shown that the order-disorder phase transitions are 
affected by the surfaces and may differ markedly from 
that in the bulk^ii^ii^i^i^. Indeed, the atoms at the surface 
layer have fewer neighboring atoms in comparison with 
the bulk, and it is well know that the magnetic order is 
strongly modified due to the lower coordination number. 
Moreover, the changes of the lattice constant near the 
magnetic surface also modify the exchange interaction 
constants. Therefore, the reduced translational symme- 
try may decrease or increase the magnetic stability. 

Depending on the value of exchange coupling constant, 
two different behaviors at the surface can be obtained: (i) 
If the coupling constant is smaller than a critical valve, 
the bulk and surface ordering occur at the same temper- 
ature. This case is called ordinary transition. Experi- 
ments at the bulk Curie temperature have shown that 
in general case, the magnetization of thin films and sur- 
faces is lower than the bulk one, and as the thickness 
of the films decreases, the magnetization also decreases. 
Ni and Fe correspond to this kind of systems^i^. (ii) If 
the exchange coupling constant is larger than a critical 
value, which is itself larger than the bulk coupling con- 
stant, the Curie temperature of thin films and surfaces 
is higher than the bulk Curie temperature. This case 
is called surface transition where the surface layer alone 
orders, while the bulk remains disordered. This remark- 
able phenomenon was observed for various compounds, 
such as for example Gd^, Tbii and NiQi^.is^ yIow- 
ever, for Gd it should be mentioned that there are also 
newer studies claiming that no separate transition at the 



surface exists (e.g. see RefsJ^ii^). Challenges in experi- 
mental studies of surface magnetism stem from the diffi- 
culties of preparing and characterizing clean, structurally 
and magnetically ordered Gd surfaces. 

From theoretical point of view, the spin-dependent 
band structure of films and surfaces have been investi- 
gated in local-moment systems in which the magnetic 
moment originates from a partially filled shell of the 
atoms and being strictly localized at each magnetic ion 
site. Kurz et ali^ studied the magnetism and the elec- 
tronic structure of Gd(OOOl) surface on the basis of den- 
sity functional theory. Schiller et ali^t^ investigated the 
temperature-dependent band structure of ferromagnetic 
semiconductor films and surfaces in the s-f model^^ . Us- 
ing a moment-conserving decoupling procedure for suit- 
able defined Green's functions, the layer-dependent mag- 
netizations and the LDOS were studied at different tem- 
peratures. 

The aim of this work is to present the single-site coher- 
ent potential approximation (CPA) for the s-f to inves- 
tigate the spin-polarized surface states in local-moment 
systems. To our knowledge, this is the first detailed study 
of the surface effects on electron-spin polarization using 
the CPA. In this work, based on the single-site CPA for 
the s-f model, we intend to develop a theory for elec- 
tronic states of a magnetic semi-infinite crystal, which 
is applicable in wide range of temperature and the s-f 
exchange coupling strength. 

The CPA is available when the scattering by the s-f 
exchange interaction is equivalent to the scattering by a 
completely random short-ranged potential. In this ap- 
proximation the multiple scattering on the single site is 
taken into consideration. The single-site CPA gives rea- 
sonable results for both weak and strong exchange inter- 
action limits, and hence is fairly good for any values of 
bandwidth, scattering potential and temperature^-. The 
single-site CPA, however, does not consider the exchange 
scattering caused by the /-spin correlation between dif- 
ferent sites, which plays an important role around Curie 
temperature in local moment systems. In fact, the collec- 
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tive mode, correlation, and/or clustering effect of local- 
ized spins are completely beyond the scope of the single- 
site CPA. 

The crystal is assumed to consist of completely filled 
and well ordered ferromagnetic layers and its surface is 
flat. Thus, the effects due to vacancies, islands, steps or 
intermixings on the magnetic properties are not consid- 
ered here. However, the effect of magnetic anisotropy will 
be taken into account in the values of coupling constant 
between atoms on the surface of crystal. 

The article is organized as follows: In section II, by ap- 
plying the single-site approximation to the s-f model, we 
present a theory to study the surface effects on the spin- 
dependent electronic states in a magnetic semi-infinite 
crystal. In appendix A, we obtain the Green's func- 
tion for a nonmagnetic semi-infinite crystal and in ap- 
pendix B, a procedure is given for calculation of the 
layer-dependent magnetization in the frame work of the 
mean-field approximation . The results of numerical cal- 
culations of the self-consistent equations for the LDOS 
and the electron-spin polarization are discussed in sec- 
tion III for two cases: ordinary and surface transitions. 
In section IV, we conclude the article with a summary. 

II. MODEL AND FORMALISM 

We consider a magnetic semi-infinite crystal which can 
be described by a single-orbital tight-binding Hamilto- 
nian with nearest-neighbor hopping t on a simple cubic 
(s.c.) crystalline lattice with lattice constant a. The 
structure is obtained by stacking layers parallel to the free 
surface and we choose the (001) axis of the s.c. structure 
to be normal to these layers. This direction is called z- 
direction hereafter and the system occupies the half-space 
z > 0. We use the s-f model which is commonly consid- 
ered as realistic for local-moment structures^^iSii^. In 
this model the following Hamiltonian is used to describe 
the magnetic semi-infinite crystal: 

Ht^H, + Hf + Hsf , (1) 
Hs=t J2 E \r,n;<j){r',n';a\ , (2) 

n,n' — 1 r,r' ,(T 

+ 00 

= — ^ ^ ^ ^ Jrn.r'n'^r.n ' ^r' ,n' i (3) 
n,n' — 1 r,r' 

+ 00 

Hsf ^-I^Yl k,";cr)(o-<T<T' • Sr,„)(r,n;CT'| . (4) 

n— 1 r,{7,a' 

Here, Hs describes the conduction electrons as s- 
electrons and gives the electron transfer energy in which 
|r, n; a) is an orbital with spin a (=1, |) at site r in the 



x-y plane and layer n (= 1, 2, • • • , -|-oo). Each lattice 
point of the semi-infinite crystal is occupied by a local- 
ized magnetic moment, represented by a spin operator 
Sr.„. The direct exchange coupling between these lo- 
calized moments is expressed by the Heisenberg Hamil- 
tonian Hf where Jrn,r'n' is the exchange coupling con- 
stant. The problem with the simple Heisenberg model 
in the form ^ is that the atoms of the topmost sur- 
face layer have fewer neighboring atoms than those of 
the bulk; thus, the magnetic order is strongly modified 
by the coordination number. This point is considered in 
calculation of the layer-dependent magnetization. H^f is 
the s-f exchange interaction between the s-electron and 
the localized /-spins where a^a' is the component of the 
Pauli matrix and I is the s-f exchange coupling constant. 
In this work, we regard a /-spin as a classical spin since 
the magnitude of the localized spin on a magnetic ion 
(S=7/2) is pretty large while the exchange interaction 
strength IS is kept finite. The summations in the above 
equations extend over the sites of the semi-infinite crys- 
tal. 

The single-electron Hamiltonian is defined as 

H = H,ff + V , (5) 

where the effective Hamiltonian H^f / which describes the 
effective medium is expressed as 

Heff ^ Hs + |r,n;fj)Bn^(r,n;f7| , (6) 

n— 1 r,fT 

and the perturbation term V is written as 

-\-oc 

V = Hsf ~^^\r,n\a)Ka{i^,n]a\ (7) 

n— 1 r,cj 

+ 00 

= ^{Y^ I"*' ^'^ ^) [-I{(^<y<y' ■ Sr,n) 

n— 1 r (7,(j' 

r, n; cr'|} (8) 

+00 

= EE^r.n . (9) 

n=l r 

Here fi„o- is the layer- and spin-dependent coherent po- 
tential and Wr,n is an isolated potential in site r of the 
ri-th effective layer. 

As in Ref.— , we apply the condition that the average 
scattering of the s-electron by the single /-spin embedded 
in the effective medium is zero. Thus we define the single- 
site t-matrix of the s-f exchange interaction as 

tv,7i^V,.,n[l- PVy^^nT'^ (10) 

where P{u}) — G^{uj — 8na) is the effective Green's func- 
tion and G°(a;) is the Green's function of the nonmag- 
netic semi-infinite crystal. Here, tr,n is the complete scat- 
tering associated with the isolated potential Wr,n- 
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Within the single-site CPA, the condition (tr,n)t/i = 
for any r in each layer, leads to the spin-dependent effec- 
tive medium where an electron is subjected to a coherent 
potential, 6„| or fi„|, according to the electron spin ori- 
entation. The spin-dependent coherent potentials are en- 
ergy (ti;)-dependent complex potentials and in this study 
described by the following equations^: 



where 

4 — 



Fr,lPS{S + 1) - ISjl + f„xfi„ J(A„|/5„t) 
1 + E^iiKl -I)- FnlIS(Ar,T/Bn^) 



(11) 



the bulk Hamiltonian. We find that the localized energies 
within the surface layer (n = 1), or at most, in the second 
layer (n = 2), i.e. the first interior layer, differ from that 
in the bulk. Here, we approximate the layer-dependent 
coherent potential, fi„o-, by the layer-independent bulk 
one, 6^, for all layers (n > 1). The difference in the 
coherent potentials is treated as a local perturbation in 
the semi-infinite Hamiltonian and applies only for the 
layers n = 1 and n = 2. Therefore, we approximate the 
effective Hamiltonian (jS]) by the following simple form: 



Fn^PSjS + 1) + 7^(1 + F„Tfi„T)(Au/g»i) 

1 + i^„T(6nT Fn^IS(AnjBnl) 

(12) 

SUS 

(1 -^ncr^icr)(l Fn—cr^n—a) ■F'ncT^n—a^^na / Hi 

(13) 

1 

(14) 



Una — Zcr^{Sj^ Z^j) B^j 



(15) 



(16) 



(17) 



In above equations (• • means the thermal average, 
Z| = -fl, Z| = — 1, and 

F^, = F,,,{lo) = ±^GL(k||;c^-fi„,) 



(18) 



where A'^n is the number of sites per layer, k|| is a wave 
vector parallel to the layers, D^^' (e) is the density of state 
for a square lattice with only nearest-neighbor hopping, 
and G'5^„(k||, is a matrix element of the Green's func- 
tion of a semi-infinite crystal which is given in Appendix 
A. 

It is clear that the Eqs. 1131)- (d?]) include only 
as an /-spin operator. Therefore, the thermal average 
for fluctuating /-spins can be easily calculated using the 
molecular-field theory. On the other hand, the existence 
of surface modifies the magnetic order in semi-infinite 
crystals. Thus, the thermal average of the /-spin op- 
erator depends on the layer index n. In appendix B we 
have derived and depicted the normalized magnetizations 
{S^)/S which depend on layer index. 

However, with increasing n from the surface layer, the 
eigenvalues of semi-infinite Hamiltonian rapidly approach 



Heff = H^ff 



+CSO 



U 



(19) 



n— 1 r,cr 



with 



"^ncr — ^na B^ 



(21) 



(22) 



Here, Una is a local perturbation for electrons with spin 
a in the rt-th layer of the semi-infinite crystal. 

The layer- and spin-dependent Green's function Gnma 
which corresponds to -ffe// is then given by 

G'„m<T(k||; w) = G"„j„(k||;w) 



where 

Glmcr 



(23) 



[1 - G?i,7.i.][l - &^,,U2„] - Gl^^u^^Gl^^u^, ' 

(24) 



[1 - G;i.^i.][1 - GO2.W2.] - Gl^^u^^Gl^^u^. ' 

(25) 



and G^TOcrC^ii;!^) = G°„(k||;ti) — fi^) is the correspond- 
ing Green's function of -ffg^j and can be obtained from 
appendix A. Therefore, using Eqs. l|^ - (P5)) the diago- 
nal elements of the Green's function for the semi-infinite 
effective crystal can be resulted from the relation 



G 



2rna 



Fr,. 



Fnaito) 



deD^^\e)Gr,na{e;to) 



(26) 



However, before determining Fna-, we should obtain fi^. 
In Eqs. (HHl-lini), by changing fi„<, ^ fi^;, -> 
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A V 



Va, Una ^ U^, Wna ^ W^, and 



-E 



1 



- fi^(c.) - e(k) 



lyjl -9(- 



(27) 



we can obtain the bulk self-energy^. Here is the 
diagonal element of the bulk effective Green's function, 
'N is the number of cells in the infinite crystal, and e(k) = 
2i(cos k^a + cos kya + cos fc^a). 

When is determined, we can obtain 6„ct using Eqs. 
(fTT]l - (|17p . However, in these equations we use of F„cr 
(Eq. (I2S1)) instead of F^a (Eq. HH])). This procedure is 
repeated until the calculation converges. After calculat- 
ing the layer- and spin-dependent Green's function Fna^ 
the LDOS in layer n can be determined by the relation 



F)na{^) = ImF„cr(ti^ + iS) 



(28) 



which should satisfy the following equation in all of the 
present numerical calculations 




■1.5 -1.0 -0.5 0.0 0.5 1.0 -1.5 -1.0 -0.5 0.0 0.5 1.0 1.5 

ra/W (i)/W 

FIG. 1: The bulk DOS as a function of energy in various 
temperatures. The dashed (solid) curve is the DOS for spin- 
up (spin-down) electrons. 



(29) 



In Eq. S is a positive infinitesimal. 

In order to study the electron-spin polarization, we 
assume that N^/Ni is equal to £)„^(tj)/D„|((jj) where 
TV-f (7V|) is the number of electrons with spin-up (down). 
Thus the degree of electron-spin polarization at Fermi 
energy lo^ and in layer n can be given by 



■P«T(^i. ) - Au(^^) 



(30) 



III. RESULTS AND DISCUSSION 

Using the approach described in the previous section, 
we have studied the spin-polarized surface states for 
IS/W— 0.35 which corresponds to the intermediate s- 
f exchange interaction. First we investigate the spin- 
dependent electronic states in different positions. 

In Fig. 1, we have shown the DOS as a function of 
energy for a s.c. crystal at T^O, 0.5, 0.8 and T > T^. 
It is clear that the largest spin-splitting between both 
spin-up and spin-down bands occurs at zero tempera- 
ture. At this temperature, the spin-up DOS is similar to 
the DOS of a nonmagnetic crystal. Since at zero temper- 
ature (completely ferromagnetic case) all the local mo- 
ments are aligned at z-direction, the magnetic system is 
periodic. As a result, the spin-up states are only shifted 
to the low energy side with no damping, and diminish at 
high energies where the coupling of the states are mainly 
antiparallel. On the other hand, the spin-down electrons 



can be flipped; thus, the spin-down band is shifted to the 
high energy region, while the bottom of the spin-down 
band extends down. With increasing the temperature, 
the local moments fluctuate and the spin-splitting be- 
tween both spin bands decreases. At T > T^, the spon- 
taneous magnetization is zero and hence the DOS are 
spin-independent . 

In Figs. 2 and 3, we have depicted the surface LDOS 
for spin-up and spin-down electrons at different temper- 
atures. The curves in Fig. 2 correspond to the case 
in which the surface orders at a temperature Tq (sur- 
face transition) larger than the bulk transition temper- 
ature Tp. In this regard, we have chosen J||=1.5 J and 
Jj^=1.15 J. The curves in Fig. 3, however correspond 
to the case in which both the bulk and the surface order 
at identical temperature (ordinary transition). In such 
case we have assumed that J||=0.8 J and J_l=0.95 J. If 
we consider the layer-dependent magnetization in Fig. 8, 
we will see that, at low temperatures, the LDOS in the 
case of surface transition is similar to those obtained in 
the case of ordinary transition. At zero temperature, the 
spin-up LDOS is similar to the surface LDOS in a non- 
magnetic semi- infinite crystaP'^. With increasing tem- 
perature from zero to T ~0.5 T^, the LDOS does not 
change considerably in both cases. The spin-splitting in 
the LDOS is a measure of the electron-spin polarization 
which depends on the layer index. Further increase of 
temperature, decreases the spin-splitting of the LDOS. 
Fig. 2 shows that, for the surface transition the spin- 
splitting is not zero at T = Tq and even at T —1.1 T^, 
however for the ordinary transition the spin-splitting will 
be zero at T = as shown in Fig. 3. 
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FIG. 2: The LDOS for the surface layer (n=l) as a function 
of energy in various temperatures. The dashed (solid) 
curve is the LDOS for spin-up (spin-down) electrons. Here, 
J|l=1.50 J and J±=1.15 J. 



FIG. 4: The LDOS for the first interior layer (n=2) as a 
function of energy in various temperatures. The dashed 
(solid) curve is the LDOS for spin-up (spin-down) electrons. 
Here, J||=1.50 J and Jx=1.15 J. 
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FIG. 3: The LDOS for the surface layer (n=l) as a function 
of energy in various temperatures. The dashed (solid) curve is 
the LDOS for spin-up (spin-down) electrons. Here, Jy =0.80 
J and Jx=0.95 J. 



FIG. 5: The LDOS for the first interior layer (n=2) as a func- 
tion of energy in various temperatures. The dashed (solid) 
curve is the LDOS for spin-up (spin-down) electrons. Here, 
J|l=0.80 J and Jx=0.95 J. 



6 



1.0 



0.8 



c 


N 



0.6 



(£0.4 
c 

^0.2 



0.0 



0.8 



g 

■to 0.6 

N 

(£0.4 

c 
a 

CO 



0.2 



0.0 



1 1 1 1 1 1 1 1 1 

(a) n=1 \\^ 




— av/W=0 
— ^av/W=-1 

\ 


\ 

1 


1 . 1 . 1 . 1 . 


1. 


1 1 1 1 1 1 1 1 i 

X 




(b) n=1 Y 




— a)p/W=0 V\ 




1 , 1 , 1 , 1 . "-i 


La 1 



0.0 0.2 0.4 0.6 0.8 1.0 

11 V 



1.2 




0.0 0.2 0.4 0.6 0.8 1.0 



FIG. 6: The degree of spin polarization for the surface layer 
(n—l) as a function of normalized temperature T/Tq: (a) 
J||=1.50 J and Jx=1.15 J, (b) J||=0.80 J and Jx=0.95 J. 
The dashed curve is the bulk electron-spin polarization. 



FIG. 7: The degree of spin polarization for the first interior 
layer (n=2) as a function of normalized temperature T/T^: 
(a) J||=1.50 J and Jx = 1.15 J, (b) J||=0.80 J and ,/x=0.95 
J. The dashed curve is the bulk electron-spin polarization. 



In Figs. 4 and 5 we have depicted the temperature 
and spin dependence of LDOS in the second layer at the 
surface and ordinary transition cases, respectively. The 
electronic states slightly resemble the DOS function for 
the infinite crystal. In this layer the spin-down band, as 
the previous electronic states, has a tail which reaches 
the edge of the spin-up band even at low temperatures. 
As the temperature increases, the spin-up (-down) LDOS 
shifts to high (low) energy side; so that, at paramagnetic 
region the LDOS for both spin-up and spin-down bands 
completely coincide to each other. It is clear that in all 
cases, there is a cusp in the LDOS which is due to the 
s-f exchange coupling strength. 

If we compare the Figs. 2(f) and 4(f) (at T=l.l T^), 
we will find that the value of spin-splitting at the surface 
layer is slightly more than the first interior layer (second 



layer). Therefore, as one proceeds into the crystal, the 
value of spin polarization begins to decrease at temper- 
atures higher than Tq. When we are far enough from 
the surface, the surface effect on the spin polarization 
completely disappears and at T > the electron-spin 
polarization becomes zero. 

In order to gain further insight into the surface ef- 
fects on the LDOS, we have shown in Figs. 6 and 7, 
the temperature dependence of the electron-spin polar- 
ization in both surface and ordinary transition cases for 
the first and second layers, respectively. In most of the 
magnetic phenomena at surfaces, interfaces and in thin 
films, such as spin-polarized electron emission and spin- 
polarized transport through a magnetic film, the elec- 
trons at the Fermi surface play an essential role. There- 
fore, the electron-spin polarization depends on the posi- 
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tion of Fermi energy. Here, in both Fig. 6 and 7, we 
have calculated the spin polarization at uJj,/W = — 1 
(near the band edge) and at oj^/W = (near the band 
center). For comparison, we have shown the electron- 
spin polarization for the infinite crystal. It is clear that, 
in surface transition case when ujp/W = —1, the spin 
polarization on the surface {n=l) and in the first interior 
layer {n—2) is higher than the bulk one. The spin polar- 
ization on the surface layer differs with the bulk one at all 
temperatures. This difference near reaches its max- 
imum value. However, due to the reduction of the spin 
polarization in the first interior layer in comparison with 
the surface layer, the difference in the spin polarization 
of this layer and the bulk appears above 0.6 T^. 

The difference of the spin polarizations between both 
layers and the bulk is originated from the strong exchange 
couphng with magnetic ions. At uj^/W = 0, the spin po- 
larizations have low values at all temperatures, because 
the difference between two spin bands is small. If we 
consider the ordinary transition case, we will see that, 
near Tq that the surface magnetization has the most dif- 
ference with the bulk one, the bulk spin polarization is 
more than the value of spin polarization at the layers 
n = 1 and n = 2. 

We applied the present model for s.c. local moment 
systems such as Eu chalcogenides. This model is also 
applicable to the case of hep Gd. However, for an atom 
at the Gd(OOOl) surface, we have six nearest-neighbor 
atoms at the surface, and three nearest-neighbor atoms 
in the second layer. Applying the molecular-field theory 
to the case of a hcp(OOOl) surface, we find that, for an 
exchange coupling ratio that satisfies J" > | J, there is 
a surface transitions^. Here, is the exchange coupling 
constant at the surface. 

It should be noted that, real surfaces are usually not 
perfectly smooth and display some degree of roughness 
due to the presence of surface defects, such as for exam- 
ple steps, islands, or vacancies. Impurities are also often 
encountered at crystalline surfaces and may be viewed 
as the source of some disorder at the surface. All these 
defects have some impact on magnetic surface quantities. 
Based on the single-site CPA, Hong'^^ studied the LDOS 
near a rough surface. The numerical results show that, 
the LDOS depend sensitively on surface roughness. On 
the other hand, Zhao et al.— using Ising-model Monte 
Carlo simulations have shown that surface roughness 
strongly affects surface magnetization, not only chang- 
ing the surface magnetic ordering temperature but also 
modifying the magnitude of magnetization in a complex 
fashion. Therefore, both the LDOS and the surface mag- 
netization are sensitive to the surface defects and con- 
sequently the resulting spin polarization is strongly in- 
fluenced by the vacancies, islands, steps and the other 
surface defects. 

The above results show that the exchange coupling be- 
tween the surface magnetic ions plays an important role 
in the degree of spin polarization for the emitted elec- 
trons from the free surface or the transmitted electrons 



through interfaces. Hence, in designing the spin elec- 
tronic or "spintronic" devices the information about the 
surface effects on the degree of spin-polarization plays an 
essential role in generation of highly polarized electrons. 



IV. CONCLUSIONS 

In summary, we have presented a new formalism to 
study the spin-polarized surface states in the local mo- 
ment systems. Based on the single-site CPA for the s-f 
model and the mean-field theory for the layer-dependent 
magnetization we have investigated how the surface influ- 
ences the spin-dependent electronic states. The effect of 
the semi-infinity of the crystal was taken into account by 
a localized perturbation method. Assuming a s.c. band 
structure for the s electrons we obtained the numerical 
results for the LDOS and the electron-spin polarization 
for various temperatures and in two different behaviors 
at the surface; ordinary and surface transitions. These 
transitions depend on the values of the exchange cou- 
pling constants between atoms on the surface and be- 
tween atoms in the layers n = 1 and n ^ 2. The re- 
sults showed that in the surface transition case, in which 
the coupling constants are higher than the correspond- 
ing in the bulk, the surface Curie temperature lies 
well above the bulk Curie temperature. Therefore, at 
<T <Tl. the electron -spin polarization on the sur- 
face has a finite value, while the bulk spin polarization 
is zero. The existence of the surface magnetization and 
hence, the spin-polarized surface states above the bulk 
Curie temperature can be useful for development of new 
electron-polarized devices. 
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APPENDIX A: GREEN'S FUNCTION OF A 
SEMI-INFINITE CRYSTAL 

The existence of surface affects the electronic eigen- 
function; therefore, due to the surface effect, the Green's 
function of a semi-infinite crystal is essentially layer- 
dependent. In this regard, for the eigenfunction of a 
semi-infinite crystal we use of a simple form as: 

|k||,fc^) = -^^^„(fc^)|k||,n) , (Al) 

where the layer-dependent coefficient Lpniki_) is given by 
the standing wave 

^n{k±) = V2sm{k±na) , (A2) 
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and satisfies the following orthogonality conditions 

= , (A3) 



(A4) 



Using Eqs. (|Aip and (|A2p . we obtain the layer- 
dependent Green's function for a semi-infinite crystal 
with a (001) surface as follows; 



G°„(k||;a;) = -^5: 



2t cos{k±a) 



(A5) 



Here N± is the number of layers parallel to the surface 
and ek|| = 2i(cos k^a + cos kyu). 

By transforming the sum over k± in Eq. (lASP to an 
integral over k±a — 0, the Green's function in a mixed 
Bloch-Wannier representation for a semi-infinite crystal 
is then given as 



-3i 



where W — 6t is the half-bandwidth and 



[/(r?)l"-™l -/(r;)l"+™l] , 

(A6) 



(A7) 



(A8) 



APPENDIX B: LAYER-DEPENDENT 
MAGNETIZATION 

Here, we like to present a generalized molecular-field 
theory to calculate the layer-dependent magnetization for 
a magnetic semi-infinite crystal. Due to a wide range 
of physical effects near surfaces, surface properties may 
be significantly different from bulk propertie o^^'^^i^^'^^ . 
These are partly taken into account by a deviation of 
exchange coupling constant Jm.r'n' in nearest neighbor 
approximation and in the neighborhood of surface. For 
simplicity, we assume in all calculations that only the 
coupling constants of the local magnetic moments at the 
surface, Jri,r'i = ^||, and also between atoms in the sur- 
face and second layer, Jri,r2 = J^_, differ from the bulk 
exchange coupling constant Jrn,r'n' — J- 

The molecular- field equations given by Cottam^S for 
spontaneous magnetizations of the semi-infinite crystals 
with s.c. structure and (001) orientation of the surface, 
are of the form 



s 



, -^11 (S!) , J± (5| 



2t{S + 1) \ J S J S 



c 1.0 



■I 0.8 
c 

O) 

2 0.6 

+-< 

c 

0) 

E 0.4 
a 

? 0.2 





1 ' 1 ' 


1 ' 1 ' 1 ' 
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FIG. 8: Layer-dependent normalized magnetization, {S^)/S, 
as a function of normalized temperature, T/T^. (a) J||=1.50 
J and Jx = 1.15 J, (b) J||=0.80 J and Jx=0.95 J. 



(^1 



5* 



S 



2t{S + 1) \ J S 



J± (Sf) , , {SI) , {SI) 



(B2) 



S 



2t{S + 1) 



{Sn 



{Sn) I {Sn 



+ 1/ 



s 



s 



s 



n > 3 , (B3) 



where Bs{x) is the Brillouin function and r 



T/T^ 



which is the bulk Curie temperature of the sample. 
Depending on the values of J\\ and J±, two different be- 
haviors at the surface can be obtained: (i) If the coupling 
constants J|| and J± are larger than J, the magnetization 
at the surface and even in the layers near to the surface 
may have nonzero values at T > Tq. Therefore, these 
layers are in the ferromagnetic phase for < T < 
{Tq is the surface Curie temperature), while the bulk 
(Bl) is in paramagnetic state, (ii) If Jii and J± are smaller 
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than J, the surface and the layers near the surface may 
have weaker magnetization than the bulk at tempera- 
tures near T^. The normalized magnetization in each 
case can be determined using the self-consistent solution 
of Eqs. (jBl[) - (|B3|) . The layer-dependent magnetization 



for both cases has been shown in Fig. 8. It is important 
to note that, the temperature dependence of the mag- 
netizations depends strongly on the boundary conditions 
imposed in the calculation. Here, the bulk boimdary con- 
ditions were imposed in 10th layer—. 
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